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Chapter 7 

Bio-Optical and Remote Sensing Observations in 

Chesapeake Bay 

Lawrence W. Harding Jr. and Andrea Magnuson 
University of Maryland Center for Environmental Science, Cambridge, Maryland 

7.1 INTRODUCTION 

The high temporal and spatial resolution of satellite ocean color observations will prove invaluable for 
monitoring the health of coastal ecosystems where physical and biological variability demands sampling 
scales beyond that possible by ship. However, ocean color remote sensing of Case 2 waters is a challenging 
undertaking due to the optical complexity of the water. The focus of this SIMBIOS support has been to 
provide in situ optical measurements from Chesapeake Bay (CB) and adjacent mid-Atlantic bight (MAB) 
waters for use in algorithm development and validation efforts to improve the satellite retrieval of 
chlorophyll (chi a) in Case 2 waters. CB provides a valuable site for validation of data from ocean color 
sensors for a number of reasons. First, the physical dimensions of the Bay (> 6,500 km 2 ) make retrievals 
from satellites with a spatial resolution of ~1 km (i.e., SeaWiFS) or less (i.e., MODIS) reasonable for most 
of the ecosystem. Second, CB is highly influenced by freshwater flow from major rivers, making it a 
classic Case 2 water body with significant concentrations of chlorophyll, particulates and chromophoric 
dissolved organic matter (CDOM) that highly impact the shape of reflectance spectra. Finally, past and 
ongoing research efforts provided an extensive data set of optical observations that support the goal of this 
project. Our SIMBIOS contribution has the following objectives: 

• To provide to SeaBASS an on-going stream of in-situ optical measurements from CB/MAB for 
use by the ocean color community in the development and validation of atmospheric correction 
and ocean bio-optical models; 

• To evaluate the performance of SeaWiFS bio-optical algorithms in CB/MAB; 

• To evaluate the performance of SeaWiFS L wv and chi a products in CB/MAB on a wide range of 
time and space scales; 

• To utilize our bio-optical database to develop a regionally-tuned algorithm for CB/MAB; 

• To produce accurate satellite chi a fields that can be used in existing primary productivity models 
(Harding et al., 2002) to improve the spatial and temporal resolution of primary productivity 
estimates for CB/MAB. 


7.2 RESEARCH ACTIVITIES 

Field Program 

Throughout our SIMBIOS support we have maintained a seasonal (i.e., spring, summer, fall) sampling 
schedule in CB proper that started in 1995 with the NSF Land-Margin Ecosystem Research (LMER) 
program Trophic Interactions in Estuarine Systems (TIES) cruises (1995-2000) and continued with the NSF 
Biocomplexity cruises (2001-2003). In addition, five separately funded Office of Naval Research (ONR) 
cruises were conducted in adjacent offshore waters including coverage to the Gulf Stream. All cruises 
collected the same suite of measurements. Vertical profiles of downwelling irradiance (Ed) and upwelling 
radiance (L u ) were preformed with a Biospherical Instruments MER 2040/2041. Discrete samples of 
particulate and dissolved fractions were collected for determination of absorption spectra of phytoplankton, 
a ph (k), non-pigmented particulates, a d (L), and colored dissolved material, a C DOM(7)- Discrete particulate 
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samples were collected for pigment analyses of chi a by fluorometry and a suite of pigments by HPLC. 
Radiometric profiles were used to calculate attenuation coefficients for downwelling irradiance (IQ) and 
upwelling radiance (K u ), and surface values of E d and Lu (E d 0" and L u 0", respectively) by simple linear 
regression. Instrument self-shading corrections were applied to each L u 0" according to Gordon and Ding 
(1992). 

In 2001 we expanded our sampling protocol to include two new radiometric instruments, the Satlantic 
MicroPro free-fall radiometer and the Satlantic Hyperspectral Tethered Spectral Radiometer Buoy (Hyper- 
TSRB). The MicroPro is designed for work in more turbid waters with a small sensor cross-section, and is 
sufficiently compact to allow deployment from small boats. The Hyper-TSRB measures surface E d and L u 
with a 120-channel detector from 400 to 800 nm. On each cruise in CB, sequential radiometric profiles 
were made with the MER (typically one cast) and MicroPro (2 to 5 casts). All casts were completed within 
~5 to 8 minutes, depending on the number of MicroPro casts that were made. The TSRB was deployed 
prior to the radiometric profiles and continuously operated while on station. Since 2001 we also have been 
using a MicroTops sun photometer from the SIMBIOS instrument pool to measure aerosol optical 
thickness. 

In 2002 and 2003 we participated in six EPA/NASA Atlantic Coast Estuaries Indicators Consortium 
(ACE INC) cruises that sampled two tributaries of CB, the Patuxent and Choptank Rivers. ACE INC 
cruises collected the same suite of measurements as the LMER-TIES and Biocomplexity cruises, with the 
exception that radiometric profiles of E d and L u were limited to the MicroPro because of vessel limitations. 

We have submitted over 3000 in-situ measurements to SeaBASS, significantly expanding bio-optical 
data for Case 2 waters (Table 7.1). Data from summer and fall cruises of 2003 will be submitted prior to the 
end of this contract. We also participated in the CHORS/Horn Point Laboratory HPLC Round Robin. We 
worked together with Jason Perl to identify and minimize methodological and instrument differences 
between laboratories such that an internally consistent pigment database could be created. 

Algorithm Evaluation/Development 

We have evaluated the performance of the empirical algorithm OC4v.4 (O’Reilly et al., 1998) and the 
semi-analytical Garver-Siegel-Maritorena (GSM01) model (Maritorena et al., 2002) for the CB/MAB, 
using bio-optical data collected during from SIMBIOS. This entailed using regional bio-optical data as 
input parameters to GSM01, specifically the chi a-normalized absorption, a ph *(k), and spectral slope of the 
combined absorption of dissolved and non-pigmented particulate matter, S C d m , to produce a regionally tuned 
version of the model that we termed GSM01-CB. We subsequently modified the native GSM01 parameters 
in the SeaWiFS Data Analysis System (SeaDAS) code to evaluate the performance of the GSM01-CB 
using satellite-derived radiances as input. 

Sea WiFS Image Analysis 

We have continued acquisition and analysis of SeaWiFS ocean color imagery from the Goddard 
Distributed Active Archive Center (DAAC) for CB/MAB. Approximately 100 clear images are available 
each year fro 1998-present. We have used in-situ data to validate SeaWiFS performance in match-up 
analyses and comparisons with continuous underway measurements. A three-year times series (1998-2000) 
has been used to determine the extent to which SeaWiFS resolves seasonal and interannual variability of 
chi a evident in the in-situ data. 


7.3 RESEARCH RESULTS 

Radiometric Comparisons 

Expansion of our radiometric measurement suite allowed us to evaluate performance of each profiling 
instrument in nearly simultaneous (i.e., within minutes) deployments. The MicroPro is more suitable for 
work in Case 2 waters than the MER due to: (1) smaller diameter (MicroPro 4.8 cm relative to the MER 21 
cm) that reduces instrument self-shading; (2) higher sampling frequency that gives greater depth resolution 
in waters where IQ is high and slight changes in the depth of the sensor have a larger effect on the 
measurements. The MicroPro exhibited clear advantages in waters with higher K d . The improved depth 
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Table 7.1: Status of the Chesapeake Bay data set submitted to SeaBASS. Numbers for 2003 are estimated 
based on cruises yet to be completed. 


Direct Measurements 


Direct 

Measurements 

Derived 

Products 

1996 

1997 

1998 

1999 

2000 

2001 

2002 

2003 

MER* 

E d (z,X) 

L u (z,X) 

E s (0 + ,X) 

Rrs( 0\X) 
K d (m,X) 
K u (m,X) 

82 

93 

83 

44 

40 

25 

39 

25 

MicroPro 

E d (z,X) 

L u (z,X) 

E s (0 + ,X) 

Rrs( 0 + ,X) 
K d (m,X) 
Ku (m,X) 

— 

— 

— 

— 

— 

49 

130 

90 

a p (0,X), a d (X) 

@ph(X) 

79 

32 

130 

66 

46 

27 

32 

25 

c * cdom(0,X ) 


79 

32 

130 

66 

46 

27 

32 

25 

Chi a (z) 


323 

709 

591 

514 

330 

189 

182 

60 

Pigments (z) 


143 

95 

104 

63 

56 

49 

32 

20 

T(z), S(z) 


— 

30 

55 

49 

43 

26 

31 

20 

MicroTops 


— 

— 

— 

— 

— 

— 

24 

15 


* MER wavebands 412, 443, 455, 490, 510, 532, 550, 650, 589, 625, 671, 683, 700 

**MicroPro wavebands 400, 412, 443, 455, 490, 510, 532, 555, 565, 590, 625, 670, 684, 700. 

E d (z,X) Downwelling irradiance depth profiles 

L u (z,X) Upwelling radiance depth profiles 

E s (0 + ,X) Incident irradiance 

a p (z,X) Total particulate absorption (300-900 nm) of discrete samples by QFT 

a d (z,X) Non-pigmented particulate absoiption (300-900 nm) of discrete samples by 

QFT with MeOH extraction 

a c dom(z,X ) Colored dissolved absoiption (190-820 nm) of discrete samples 

Chi a (z) Determination of chi a for discrete samples by Turner fluorometry 

pigments (z) HPLC analysis for a suite of pigments 

T(z), S(z) SeaBird temperature and conductivity profiles 


Derived Products 


Rrs(0\X) 

Remote sensing reflectance from calculated from radiance/irradiance 
profiles 

K d (m,X), K u (m,X) 

Attenuation coefficients for E d (z,X) and L u (z,X) for the surface mixed layer, 
calculated from radiance/irradiance profiles 

®ph(z,X) 

Pigmented particulate absoiption (300-900 nm), calculated as a p (0,X)- 
a d (0,X) 
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3 

In-situ chi a (mg m ) 

Figure 7.1: Comparisons of in-situ chi a to: (a) mean chi a estimated from OC4v.4 applied in-situ RRS; (b) 
SeaWiFS retrievals of chi a. SeaWiFS values represent the mean of a 3x3 pixel box centered on the station 
location for in-situ sampling. A chi a match-up was considered valid if the satellite observation and in-situ 
measurement occurred on the same day, and half (5 out of 9) pixels in a box centered on the latitude and 
longitude of the sampling station returned valid estimates. The original match-up dataset contained over 
1400 chi a observations. Mean ratio = SeaWiFS chi a / in-situ chi a. 
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Figure 7.2: RRS555 and chi a from SeaWiFS for a CB scene, 15 October 1999. The white lines 
superimposed on these images give the ship track for concurrent measurements that included bio-optical 
measurements. Vertical line at ~ 150 km marks the time of the satellite pass. 

resolution of the MicroPro (~25 measurements per meter) compared to the MER (~10 measurements per 
meter) compensated for the shallow optical depths in CB where the depth range available for the 
extrapolation interval was limited to the upper 4 m (or less). Measurements of E (1 0 from the MicroPro were 
much easier to reconcile with surface irradiance measurements (E s ) than were those from the MER, which 
exhibited an increasing bias as K d increased. Measurements of L u 0" from the MER and MicroPro also 
showed increased disagreement as K d increased. Based on advantages of the MicroPro over the MER in 
comparisons of E d 0 1 and E s , we feel confident that the MicroPro measurements of L u 0" are more accurate 
than the MER for the Case 2 waters in which we work. Lastly, the increased depth resolution of the 
MicroPro reduced the effects of small changes in the selection of the depth interval on the calculated K d , 
K„, E d 0", and L u 0". The MicroPro also allowed us to expand our radiometric sampling to more turbid waters 
of the tributaries, and has improved measurement capabilities in the upper reaches of CB. In addition, these 
comparisons of the MER and MicroPro indicated strong agreement in the clearer waters of the mesohaline 
and polyhaline CB, confirming confidence in our multi-year time series of radiometric measurements with 
the MER in the CB/MAB. 


88 




SIMBIOS Project Annual Report 



c/a 




& 



50 100 150 200 

Distance along transect (km) 


250 


Figure 7.3: (a) Comparison of in-situ and SeaWiFS chi a; (b) SeaWiFS RRS510 and RRS555; (c) ratio of 
SeaWiFS RRS 510/555 along a ship track in the mesohaline to polyhaline CB. 

Characterization of Bio-optical Properties in Chesapeake Bay 

Distributions of inherent optical properties in CB/MAB showed the strong influence of freshwater flow 
that complicates retrievals of chi a from satellite ocean color measurements. a p h(440), acDovi(440), ad(440) 
ranged from 0.1 to 1.0 m" 1 in CB, and were an order of magnitude lower in the MAB. Whereas a d (440) 
dominated in upper CB, the magnitudes of a ph (440), acDOM(440), a d (440) were comparable in the mid- and 
lower CB and inshore MAB. a C DOM(440) dominated water column absorption in the offshore MAB. The 
distribution a CDO M(440) was conservative in CB proper, due to its predominant terrestrial origin. a p h(440) 
and a d (440) were subject to much higher spatial variability. The distribution of dissolved and particulate 
materials regulated the magnitude and spectral shape of apparent optical properties such as the attenuation 
coefficient for downwelling irradiance, K d , and the remote-sensing reflectance, RRS. 


OC4v. 4 Algorithm Evaluation 


The distribution of inherent optical properties had a pronounced effect on the performance of OC4v.4 
in CB/MAB (Fig. 7.1a). OC4v.4 returned chi a estimates for CB that had a small positive bias and high 
uncertainty, reflecting that a ph , a CD0 M, ad did not co-vary. A larger positive bias of chi a occurred in the 
MAB where too much of the total absorption was attributed to phytoplankton. Match-ups of SeaWiFS chi a 
with in-situ data showed a larger positive bias in CB than observed for OC4v.4 applied to in-situ 
radiometric data, suggesting SeaWiFS-derived radiances introduced additional errors (Fig. 7.1b). Whereas 
the fourth SeaWiFS reprocessing lessened the frequency of negative LWN, comparisons of in-situ and 
satellite-derived LWN as part of this work and others have shown that low and negative LWN from 
SeaWiFS remain a problem in coastal and estuarine waters (Flarding et ah, 2003; O’Reilly and Yoder, 
2003). SeaWiFS LWN at ‘blue’ wavebands are low in CB, causing the return of higher-than-observed chi a 
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values by OC4v.4. The lack of increased bias in chi a for match-ups from the MAB suggested errors in 
SeaWiFS chi a could be attributed primarily to OC4v.4 rather than errors in atmospheric correction. 

Resolution of Inter annual, Seasonal and Event-scale Variability by SeaWiFS 

High uncertainties we detected in match-up analyses of in-situ and SeaWiFS chi a restricted the 
usefulness of SeaWiFS data for CB/MAB, stimulating additional analyses to lessen time-space aliasing we 
believed was responsible. High frequency and small-scale variability of biomass typical of estuaries 
virtually assures time-space aliasing of satellite observations and in-situ measurements. To minimize 
aliasing, we compared continuous shipboard measurements of chi a were to SeaWiFS chi a. along ship 
tracks to evaluate the extent to which SeaWiFS resolved variability on small spatial scales. We found 
SeaWiFS accurately captured small-scale (tens-of-kilometers) variability in the mesohaline and polyhaline 
regions of CB (Figs. 7.2 and 7.3) (see Harding et al., 2003). Increased noise in chi a with distance north 
corresponded to shaip gradients of RRS. We suspect that the magnitude of the RRS signal often approaches 
the limits of sensitivity for the SeaWiFS sensor in turbid waters with high K d typical of mesohaline and 
oligohaline CB, thereby limiting the accuracy of chi a retrievals. This gradient of RRS was persistent in the 
imagery and was suggested in field radiometric measurements. Comparisons of several transects suggested 
axial movement of this gradient of RRS could determine the areal extent of accurate chi a retrievals. 

Our evaluation of the 1998 - 2000 SeaWiFS chi a time series demonstrated that SeaWiFS accurately 
represented seasonal and interannual variability of chi a revealed by in-situ data (Fig. 7.4). Previous 
research has shown that the timing, position, and magnitude of the spring bloom are regulated by 
freshwater flow from the Susquehanna River, supplying nutrients and suspended matter to CB (Malone et 
al., 1988; Harding, 1994). We found that positive deviations of flow in 1998 resulted in a large spring 
bloom that extended to the polyhaline CB, shown in both satellite imagery and in-situ data. Measurements 
of chi a in 1999, a year of lower-than-average freshwater flow, showed a smaller bloom that was confined 
primarily to the mesohaline CB. Lastly, average freshwater flow in 2000 supported a bloom of intermediate 
magnitude that extended throughout the mesohaline and part of the polyhaline CB. As observed in transect 
data, the agreement of SeaWiFS and in-situ chi a improved with distance down the Bay. These 
comparisons with in-situ data confirm that SeaWiFS chi a retrievals are reasonable and can significantly 
augment the temporal resolution of shipboard measurements. 

The benefit of increased temporal resolution afforded by satellite observations was demonstrated in the 
three-year time series of chi a (Fig. 7.5). Mean, daily chi a for a section of polyhaline CB indicated several 
doublings over a period of days in summer that were not captured by in-situ observations. We hypothesize 
that these doublings could be associated with short-lived, but dense blooms of dinoflagellates that have 
been observed before during stratified conditions in summer. The sedimentation of organic matter from the 
spring bloom fuels at least 50% of the primary productivity in surface waters (Kemp and Boynton, 1984). 
Transient destratification of the summer pycnocline provides the nutrients to support these short-lived 
summer blooms (Malone et al., 1986, 1988). The higher frequency of these events in 1998 and 2000, 
contrasted with 1999, may indicate a connection between these transient blooms in summer and the 
magnitude of the spring bloom. Further study is needed to test this hypothesis. 

Regional Parameterization of GSM01 

Semi-analytical models such as GSM0 1 have the potential to distinguish multiple constituents of the 
water column that absorb and scatter light. This attribute is crucial for accurate retrievals of chi a in 
optically-complex waters. The original parameterization of GSM01 was based on in-situ observations from 
Case 1 waters and optimized for performance in the global ocean (Maritorena et al., 2002). We 
hypothesized that regional parameterization was necessary to improve the accuracy of chi a retrievals in 
CB/MAB. 

We evaluated our extensive data set of in-situ observations to derive regionally and seasonally (when 
appropriate) values for a ph *(k) and Scdm to produce a regional version of the model, GSM01-CB (Table 
7.2) (see Magnuson et al., 2003). Regional differences of a ph *(k) were significant across the study area 
from the upper, oligohaline CB to the offshore MAB. a ph *(k) increased with distance offshore. This trend 
of increasing a p h*(L) can be explained by two factors: decreases of pigment packaging and increases in 
concentrations of accessory pigments (Yentsch and Phinney, 1989; Bricaud et al., 1995). 
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Figure 7.4: Monthly mean chi a for the polyhaline (I, II, south of 37.8°N) and mesohaline (III, IV, 37.8 - 
38.8°N) regions of CB from in situ (closed symbols) and SeaWiFS (open symbols) data for 1998-2000. 


We also found significant seasonal differences in mean a p h*(k) in CB that we attributed primarily to 
changes in floral composition and associated changes in pigment packaging. Highest a p h*(k) occurred in 
spring when phytoplankton in CB consist primarily of large cells (i.e., centric diatoms) (Malone et al., 
1988; Marshall and Nesius, 1996). Lowest a p h*(k) occurred in summer when phytoplankton consist of a 
more diverse flora of smaller, flagellated cells (Malone et al., 1991; Marshall and Nesius, 1996). Scdm 
exhibited an increasing trend with distance offshore, reflecting the decrease in the influence of terrigenous 
CDOM (Blough et al., 1993; Nelson and Guarda, 1995; Vodecek et al., 1997). Small, but significant 
differences in S c( j m along the axis of CB reflected changes in the relative contributions of a C noM and ad to 
acdm- 

Evaluation of chi a from OC4v.4, GSM01, and GSM01-CB for CB, inshore MAB, and offshore MAB 
showed that GSM01-CB improved upon OC4v.4 in all regions (Table 7.3). 
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Day of Year 


Figure 7.5: Time-series of mean, daily SeaWiFS (open symbols) and in-situ (closed symbols) chi a for a 
section of CB from 37.3° to 37.8°N. 


92 


Number of Occurrences 


SIMBIOS Project Annual Report 


25 

20 

15 

to 

5 

0 

1000 

800 

MW 

400 

200 

0 

I 

1000 

800 

000 

400 

200 

0 

1000 

MXl 

600 

400 

200 

0 


29 .tyr • 2 May 2000 
Clwtapmk# Bay 
M,:m 17. 17" 
Median 19.855 
si) 27.558 
n - 40 


In -situ 


(a) 


^11 




LiL 


20 




O 

© 


• l i i i i l l l l l l t 

4 May 2000 
Clh&apeafc* Bav 
Mean 14 06 
Mcdbtt 13.200 
SI) 7.193 
ii 2277 


OC4v.4 


(J) 


jJ 


© 

o 


I I! IT I I ! I II I ! 1 

4 May 2000 
Clte*a peake Bay 
Mean 5 451 
Median 5.184 
SB 3.293 
n 1694 


(&) 


r i m 

< iSMO I 


c 

6 


4 May 2000 
Chesapeake Bay 
Mean - 18 861 
Median - 14.265 
SI) - 20.653 
n- 1770 


GSM01CB 


0 ) 


—fTfrfl 11 11 


May - Sep! 
Mere MAB 
Mean 

Median 3 48CI 
SI) 2.991 
n 35 

ln-iitu 

llvl . . . rrrJlflJ 1 

in 


40 


hi-situ 


JO 


10 


ok.fillllfl 


»iii 
M ay • Scpl 
OtlUwrc M AD 
Mean 0.436 
Median 0.280 
SI) - 0.280 
n- 134 



Log chi a 

Figure 7.6: Histograms of in-situ and SeaWiFS chi a observations for CB, inshore and offshore MAB for 
image from 4 May 2000: (a-c) in-situ, (d-f) SeaWiFS OC4v.4, (g-i) SeaWiFS GSM01, and (j-1) SeaWiFS 
GSM01-CB. In-situ data from the MAB were not available for this image. The inshore and offshore MAB 
histograms (b-c) reflect all field data available for those regions during May-September 1996-1998. 
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Table 7.2. Parameter values used in GSM01 and GSM01-CB. 


Region 

Season 

o-ph (412) 

a p h(443) 

a P h (490) 

a ph ’(510) 

&ph (555) 

o- P h (671) 

5>cdm 

'1 

GSM01 


0.00665 

0.05582 

0.02055 

0.01910 

0.01015 

0.01424 

0.02061 

1.03373 

GSM01-CB 

Spring 

0.02119 

0.02509 

0.01282 

0.00910 

0.00427 

0.02087 

0.01218 

0 

Upper 

Bay 

Summer 

0.02653 

0.02979 

0.01655 

0.01208 

0.00470 

0.02122 

0.01218 

0 


Fall 

0.02259 

0.02573 

0.01372 

0.01019 

0.00376 

0.02066 

0.01218 

0 


Spring 

0.02001 

0.02212 

0.01279 

0.00974 

0.00449 

0.01588 

0.01385 

0 

Mid- 

Bay 

Summer 

0.03345 

0.03900 

0.02318 

0.01664 

0.00609 

0.02285 

0.01385 

0 


Fall 

0.02758 

0.03080 

0.01826 

0.01438 

0.00691 

0.02112 

0.01385 

0 


Spring 

0.02001 

0.02212 

0.01279 

0.00974 

0.00449 

0.01588 

0.01330 

0 

Lower 

Bay 

Summer 

0.03345 

0.03900 

0.02318 

0.01664 

0.00609 

0.02285 

0.01330 

0 


Fall 

0.02758 

0.03080 

0.01826 

0.01438 

0.00691 

0.02112 

0.01330 

0 

Inshore 


0.07123 

0.08843 

0.06024 

0.04072 

0.01693 

0.03815 

0.01236 

0 

Offshore 


0.11331 

0.14678 

0.09832 

0.06048 

0.01920 

0.04349 

0.01646 

1 


Table 7.3: Comparison of the mean ratio (±SD) of modeled and measured chi a and a cdm (443), and 
estimated b bp (443) for Chesapeake Bay, inshore MAB, and offshore MAB. Mean ratio — modled/measured 



Chesapeake Bay Chi a 

Inshore MAB Chi a 

Offshore MAB Chi a 

GSMOl 

1.99 ±2.21 

1.84 ±0.92 

1.31 ±0.68 

GSM01-BAY 

1.10 ± 1.08 

0.59 ±0.57 

0.51 ±0.29 

OC4v.4 

1.20 ±0.71 

3.12 ±2.05 

2.40 ± 1.03 


Chesapeake Bay a cdm (443) 

Inshore MAB a cdm (443) 

Offshore MAB a cdm (443) 

GSMOl 

1.24 ±0.58 

1.01 ±0.54 

0.87 ±0.61 

GSMOl -BAY 

0.92 ±0.40 

1.12 ±0.42 

0.95 ±0.67 


Chesapeake Bay bb p (443) 

Inshore MAB bb p (443) 

Offshore MAB b bp (443) 

GSMOl 

0.94 ±0.65 

2.28 ± 1.35 

1.43 ±0.88 

GSMOl -BAY 

1.35 ±0.91 

2.88 ± 1.64 

1.34 ± 1.26 


The local parameterization was essential to improve performance in CB. Both GSM01 and GSM01-CB 
showed improvement over OC4v.4 in the MAB where they specifically accounted for high absorption due 
to CDOM. Analysis of SeaWiFS products using each model supported results of evaluations using in-situ 
RRS as inputs, specifically that: (1) one or both of the semi-analytical models performed as well or better 
than OC4v.4; (2) both semi-analytical models clearly improved upon OC4v.4 in the offshore MAB (Fig. 
7.6). Disagreement between the in-situ observations and model estimates for the inshore MAB likely reflect 
a sampling bias as in-situ observations were concentrated near the Chesapeake and Delaware Bay mouths 
where chi a are generally high, whereas mean chi a from SeaWiFS represented the entire inshore region. 
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Disregarding in-situ values, both semi-analytical models returned much lower mean chi a for the inshore 
MAB than did OC4v.4, and GSM01-CB chi a were lower than GSM01 chi a, consistent with results 
obtained using in-situ data as input to the models (Table 7.3). The semi-analytical models return estimates 
of a cdm (443) and bb,,(443) in addition to chi a. Estimates of a cdm (443) were closer to in-situ values than were 
estimates of chi a (Table 7.3). Meaningful evaluation of the bb p (443) product was not possible due to the 
lack of direct measurements of scattering. Relatively high uncertainties remained for all models, reflecting 
the difficulty of working in bio-optically complex waters. It is notable that the performance of GSMO 1 was 
surprisingly good, considering the parameterizations differed significantly from bio-optical properties of 
CB. The relatively strong performance of GSM01 in CB likely reflects the robust tuning procedures used 
by Maritorena et al. (2002) to derive global parameters for GSM01. 

7.4 CONCLUSIONS 

It is useful and necessary to assess the performance of satellite products at various stages to recognize 
limitations, as well as to demonstrate strengths and identify the types of applications for which current 
products are suitable. We have made use of an extensive database of bio-optical measurements to provide 
several different evaluations of SeaWiFS performance. We have identified several factors that inhibit the 
accuracy of the standard chi a product, including: (1) use of an empirical algorithm that does not account 
for all of the absorbing and scattering components of the water column; (2) atmospheric correction errors in 
coastal waters resulting in low satellite radiances; (3) RRS signals approaching the lower sensitivity limits 
of the sensor. Nonetheless, the importance of these factors lessened, and SeaWiFS recoveries improved, 
with increasing water clarity. We found that in mesohaline and polyhaline CB, SeaWiFS accurately and 
reliably represented seasonal and interannual variability of biomass, resolved event-scale phenomena that 
were missed in our field sampling, and resolved small-scale (i.e., tens-of-kilometers) regional variability. 
We are encouraged by our initial attempts to improve chi a estimates by tuning a semi-analytical model for 
regional conditions. We are exploring modifications to the parameterization that could lessen the use of 
abrupt transitions in parameter values at seasonal and regional boundaries. Improvements in SeaWiFS 
radiances also will be necessary for maximal performance of semi-analytical models in Case 2 waters. 
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